Spin Lifetimes in Quantum Dots from Noise Measurements 
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We present a method of obtaining information about spin lifetimes in quantum dots from measurements of 
electrical transport. The dot is under resonant microwave irradiation and at temperatures comparable to or 
larger than the Zeeman energy. We find that the ratio of the spin coherence times T1/T2 can be deduced from 
a measurement of current through the quantum dot as a function the applied magnetic field. We calculate the 
noise power spectrum of the dot current and show that a dip occurs at the Rabi frequency with a line width given 
by l/r, + 1/7,. 

PACS numbers: 85.35.-p, 73.63.-b, 72.25.-b, 72.70.+m 



Electron spins are promising candidates for qubits jltl. 
They constitute a natural two level system but in contrast with 
nuclear spins they can be highly polarized simply by using 
routinely available temperature and magnetic field conditions. 
Additionally one can exploit the mobility of electrons, allow- 
ing spin manipulation and measurement via charge currents, 
e.g. through tunnel barriers. 

For applications in quantum information processing it is es- 
sential that spins retain phase information for as long as pos- 
sible. This is usually quantified by the spin coherence time 
T2. Single spin Totimes have been determined optically 
and electrically Bt] . Petta et al. jstl encoded a qubit in the 
singlet and triplet states of two electron on a double quantum 
dot and demonstrated coherent manipulation, while Koppens 
et al. 10] used pulsed microwaves and top gates to demonstrate 
coherent oscillations of a single spin in a double dot. 

A method for determining spin coherence times using elec- 
trical transport through a quantum dot in the stationary state 
has been suggested by Engel and Loss jif], but this method 
can only determine coherence times up to an upper limit that 
is related to the temperature of the contacts. The study of 
current-current correlations has emerged as a tool to detect co- 
herence properties of quantum systems such as quantum dots 
101 • In nanomechanical resonators the noise power spectrum 
has been used to determine the oscillator occupation number 
and quality factor I9I. [loll . The noise power spectrum of a 
current interacting with a charge qubit (Tl, 12, ll. and 
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the current through a quantum dot under microwave radiation 
have been thoroughly analyzed 1 15 , 3 12 1 ■ 

In this Letter we suggest a method for measuring spin life- 
times using current-current correlations. At experimentally 
accessible combinations of temperature and magnetic field, 
where the Zeeman splitting is comparable with or larger than 
the thermal energy, the measurement of arbitrary intrinsic spin 
lifetimes becomes possible. To be able to detect the influ- 
ence of the finite electron spin lifetime on electronic transport 
through the dot, the dwell time of an electron on the dot has 
to be comparable to or longer than the spin relaxation times 
we wish to detect, e.g. through a suitable choice of tunnel 
barriers. We will henceforth assume this to be the case. 




Figure 1 : (Color online) The level configuration on the quantum dot 
and the distribution functions in the leads. The Fermi level of the 
electrons in the leads is the same for spin up and spin down, but in the 
rotating frame the energy of the electrons depends on their spin ori- 
entation, and the microwave frequency coq determines the difference 
in energy. Also the energy difference of the single electron states on 
the dot is reduced in the rotating frame from the Zeeman splitting A, 
to the detuning Aq, while the energy of the two electron state is sep- 
arated from the single electron states by the charging energy U. The 
bias voltage Vjj is applied symmetrically. The gate voltage Vg shifts 
the energies of the single electron levels by eVg and the two electron 
level by 2eV^. The resonant Rabi frequency fi couples spin up and 
spin down. 



Our model consists of a quantum dot coupled weakly to two 
leads. The full Hamiltonian after a rotating wave approxima- 
tion reads H = Ho + He + Hj, where Ho is the Hamiltonian 
of the dot. He the Hamiltonian for the leads and Hj describes 
the tunneling between the leads and the dot. We choose to 
model the quantum dot by including three many body states 
on the quantum dot: the single electron states for spin up and 
spin down respectively and a singlet two electron state. We as- 
sume that the triplet two-electron states are separated from the 
singlet state by a gap that is large compared to temperature and 
bias voltage. We also do not take any additional single elec- 
tron levels into account. In a magnetic field B only the single 
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electron levels are split by the Zeeman energy A- = gpLsB; the 
singlet is unaffected. Under microwave irradiation with fre- 
quency (Oq and after a rotating wave approximation the Hamil- 
tonian for the dot can thus be written 



Hn = 



(die 



where d] and ds, s =t,i, are the creation and annihilation 
operator for dot electrons in spin up and spin down state re- 
spectively and Us = d]ds. U is the mutual Coulomb repul- 
sion energy between two electrons on the dot, D. the reso- 
nant Rabi frequency, arising from an oscillating magnetic field 
in x-direction, e.g. in the electric field node of a microwave 
cavity, and Aq = A^ — hcoo is the detuning of the applied mi- 
crowave field . We have also included the action of a back 
gate that shifts the single electron energies by the amount eV^ 
(for the energy level diagram see Fig.[Tl). The Hamiltonian for 
the leads is given by 

He = ^^kas^las'^kas, 
kas 

where c[^^,Ckas, with a = L,R and s =1, |, are the operators 
for the lead electrons with momentum k. Tunneling from the 
leads to the dots is given by 

Hrit) = Y^tacla,{t)d, + h.c.. 

kas 



where are the tunneling amplitudes, which in the follow- 
ing we assume to be independent of the electron energy and 
equal for both leads. We have included the time depen- 
dence resulting from the transformation to the rotating frame 



in the lead operators, introducing c]^Q.|(r) ^ 
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To describe electronic transport through the dot we adapt a 
Markovian master equation approach llSll . The reduced den- 
sity matrix of the dot is defined as p(f) = Tr£ where 
Tr£ [. . .] denotes the trace over the environment and R is the 
the density matrix of the full system. Its evolution is given 
by the master equation dp {t)/dt = J(f ■ p (f )i where the super- 
operator J(f consists of three parts: ,J€ = J^o + J^^ + 
The part describing the free evolution is given by 

where X can be any operator. Renormalization effects due 
to interaction with the environment are included in the free 



evolution (see fllSll V We can also define the free propagator in 
Fourier space 



The influence of the electrons tunneling from the leads is de- 
scribed by the super-operator 
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Figure 2: (Color online) (a) The diagonal elements of the density 
matrix in the energy basis as a function of the detuning for differ- 
ent ratios of the decay rates Ji/Yi and Ye = lO^^Oo, y\ = IO^^coq, 
£1 = 10"'*ft)t), kgT = h(OQ, eVsd = fiCOo and Vg = 0. The dotted lines 
correspond to Yi/Yi = 50, the dashed to =10 and the solid lines 
to 72/71 = 1- (b) Diagram illustrating the influence of the different 
decoherence mechanisms on the spin dynamics. 



where = 27rg{Ef ) \ti\ , g{Ep) being the density of states 

i and Sij = — 1 other- 



at the Fermi energy, Sij = 1 for ; 

wise. ML/ff,T/i ^iL/RTf^(Oo/2, IIl/r = li±eV,a/2, g^\e) = 
gi2(e) ^ !-/(£), g^\£) = g^^{£) = /(£), /(e)being flie 

Fermi function, jz/v^'"* X = d^X and s^P'^ X = Xds, and sim- 
ilarly for the hermitian conjugate. The use of the free propa- 
gator instead of the self consistent inclusion of the full propa- 
gator does not alter the result as long as the distribution func- 
tions g'^ are smooth on the scale where the full propagator is 
peaked. In addition to tunneling, we include a phenomeno- 
logical description of the intrinsic spin decay, given by 



The first part describes an energy relaxation 
process leading towards thermal equilibrium 



■X. 
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X 



Y+ {a+G-X +Xa+a- - g+Xg-) 



71 = 1/n 



7- {(7-a+X + Xa-(7+ - a-Xa+) with 7+ + j- 
and Y+IY- = cxp(A,./A:Br) where C7+ = d^d 

The second part is a pure dephasing pro 



: d\d\^ - 

A2) 



- d^d^f. 



a = djdi and 



cess Jti^^' • X = (72 - 71 /2) (X - a-Xa-) , where 72 = 1 /T2. 
We will be interested in the electronic transport through the 
dot in the stationary state. The stationary state density matrix 
can be derived from ■ p^t ~ 0. This equation was solved 
analytically using Mathematica. For small detunings 
|Ao| ^ hcoo and hcoo < kgT the density matrix elements in 
the energy basis have the following approximate dependence 
on the detuning 



n (K \ li\n I A P'^^ ^ + Pu'^Yi/Yi 
P,7(Aoj = (;|p.vr \l) = .7 , ^T.. , (1) 



A2+i22y2/7i 



where p,^~' denotes the density matrix element far away from 
resonance and p\f^ the density matrix element on resonance. 
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Figure 3; (Color online) (a) The current / as a function of the detun- 
ing relative to the current at Aq = without microwave irradiation 
for different gate voltages and eV^^ = hcoo, kgT = hcoo, Q. = 10^'*CUo, 
Ye = lO^^CBo, 71 = 72 = 10^^ (Bf). When the magnetic field matches 
the microwave frequency the current shows a peak or a dip depend- 
ing on the applied gate voltage, (b) The current at resonance, Iq, 
relative to the current without microwave irradiation and the current 
as a function of the gate voltage, otherwise same values as in (a). 




Figure 4: (Color online) The peaked part of noise power spectrum, 
5s(ft)), relative to the noise floor around the Rabi frequency at reso- 
nance (A,. — hW()) and for different rates Yi and Vg = , eV^d = Ticoq, 
kgT — h(0(), = 10"'*coo, Ye = lO^^COo and a fixed ratio Yi/Yi = 10. 
The solid line corresponds to 71 = 5 ■ lO^^ftlo, the dashed line to 
7i = lO^^GJo and the dotted line to 71 = 2 ■ IO^^gjq- 



The dependence of the density matrix elements on the mag- 
netic field for Vg = is shown in Fig.|2ja). Far off resonance 
the ratio between the spin up and spin down populations is 
thermal while on resonance the three populations equalize. 
The width of the transition from the thermal state to the equal- 
ized state is given by 72/71. This can be understood with 
the help of Fig. |2b); The spin precesses around the direction 
of the effective magnetic field, B/gjis = (ii,0,Ao). For the 
parts of the precession closer to the z-axis the 71 processes 
drive the spin towards a thermal state along the z-axis of the 
Bloch sphere. Further away from the z-axis the 72 processes 
dominate, tending to reduce the x and y component to zero. 
The precession mixes both processes and the relative rates of 
both relaxation processes determines the stationary state. For 
72 > 7i the spin reacts to a different direction in the magnetic 
field further from resonance, resulting in a wider transition 
from thermal to equalized population. 

The current through the dot in the stationary state can be 
expressed as / = eTio ■ Pst\ , where Tio [. . .] denotes the 
trace over the dot degrees of freedom, with the current super- 
operator given by 



J as 



,i¥'jJJ^ 1,2. 



For equal tunneling rates for both leads it is sufficient to con- 
sider only the current through one of the tunneling contacts 
(L or R). In the region of small detunings, |Ao| <C Vsd, T, the 
current can be written approximately as / = /nPii + hiPii + 
'33P33, where the currents ii^j^,k^ 1,2,3 depend on V^d, Vg, T 
but only weakly on the detuning. The magnetic field depen- 
dence of the stationary current for different gate voltages is 
shown in Fig. [3] Close to the resonance the current shows a 
peak or a dip, depending on the bias voltage. The dependence 
of the current on the detuning follows from the detuning de- 



pendence of the density matrix elements as 

A2 + i22^/7i ' 

where /„o is the current far away from the resonance and Iq 
the current at resonance. It is therefore possible to extract the 
ratio of the dephasing rates 71 /72 from the width of the current 
peak if the Rabi frequency D. is known. 

The noise power spectrum consists of a constant noise floor, 
Sf, and two features, Ss{a)), sitting on top of that noise floor, 
a peak at zero frequency and a dip at the Rabi frequency. A 
similar noise power spectrum was discussed by Aguado and 
Brandes for a charge qubit in II12II . The floor is given by 



where we introduced the diffusion super-operator 
& = Ye fdeY,\g'\£-^asW^i'^-%{e)-K 

J as 

+ gi'{e + iias)^^} ■ %{e) ■ ^if ] , / ^ y, i,j = 1,2. 

and the features in the noise power spectrum are described by 
the term 

with 8 ^ = J' —I and the evolution super operator 



This frequency dependence of the features in the noise power 
spectrum can be evaluated at resonance to give in the positive 
half of the spectrum 



Ssi(o) = So 



Sn- 
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where Sq jq are the peak heights at zero and at the Rabi fre- 
quency. The peak and the dip have both a Lorentzian Hne- 
shape. The peak at zero frequency stems from elastic pro- 
cesses on the dot, its width is determined solely by the tun- 
neling rate jg, while the dip appearing at the Rabi frequency 
is due to processes involving energy exchange on the dot and 
its width depends on the average of the two spin decoherence 
rates, Jq, ^ (71 +72)/2. 

We have now all the ingredients to deduce the spin life- 
times T\ and 72. We can envisage the following measure- 
ment procedure: A scan of the gate voltage without an ap- 
plied microwave field will produce a series of resonances in 
the current. Selecting one of the resonances and scanning the 
magnetic field while applying microwave radiation will give 
information about the presence of an active spin on the dot 
and the width of the resulting peak gives the ratio T\/T2- Fi- 
nally adjusting the magnetic field to be exactly on resonance 
and measuring the noise power spectrum around the Rabi fre- 
quency we can obtain \/T\ + 1/72, thus making it possible to 
deduce T{ and T2. 

Why is is necessary to measure the current noise spectrum 
at finite frequency? Engel and Loss showed that measuring 
only the current through a quantum dot as a function of the 
gate voltage can already give information about T2, but to re- 
solve a certain line-width in the dot the distribution functions 
of the leads have to be sharper than the line-width of the level 
under consideration. The typical sharpness of the distribu- 
tion function is given by temperature. The achievable tem- 
peratures in today's dilution refrigerators of ~ lOmK corre- 
spond to a resolvable line-width of ^ 100 MHz, while values 
for Ti j2 in group IV semiconductors correspond to line-widths 
of < 100 Hz 1191]. Our scheme will work even at temperatures 
high compared to the intrinsic line-widths of the dot levels 
and works well at temperatures comparable to or larger than 
the Zeeman splitting of the single electron levels, although for 
temperatures much larger than the Zeeman splitting signal to 
noise ratio deteriorates. The method operates in the steady 
state and no pulsed gates are necessary. 

What are the limitations of our method? To detect a cer- 
tain intrinsic decay rate the tunneling rate has to be smaller 
than that rate. For expected decay rates of the order of 100 Hz 
that would correspond to a current of 0.01 fA. Currents of 
this magnitude are measurable using electron counting tech- 
niques 1I20I1 . At the same time we need to detect the noise 
power spectrum around the Rabi frequency, typically of the 
order of MHz. To resolve the peak the charge counting has 
thus to work at rates of 10 MHz or higher Radio frequency 
(RF) single electron transistors or RF quantum point contacts 
have been shown to work as charge counters up to 1 GHz and 
should make such measurements possible. 

In our modeling we phenomenologically introduced the 
spin decay rates as constants. A more detailed analysis of 
the decoherence mechanisms shows that the coherence time 
72 will, in principle, depend on the Rabi frequency. Assum- 



ing a Markovian bath the decay rate is directly related to the 
spectral density of the bath at the Rabi frequency. Changing 
the Rabi frequency by altering the microwave intensity should 
make it therefore possible to probe the spectral density of the 
bath. A dependence of the coherence time on microwave in- 
tensity has been measured in 10]. With this simple model we 
hope to have shown the usefulness of noise measurements in 
determining spin lifetimes. 
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